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INTRODUCTION 


There  ere  several  competing  factors  that  limit  the  performance  of  the 
monopulse  azimuth  estimation  accuracy  in  a  noise  and  interference  background. 
The  first  is  due  to  the  fading  phenomenon  that  results  in  a  degradation  in  the 
signal-to -noise  ratio  (SNR)  as  the  target  and  interference  signals  approach 
the  out  of  phase  condition.  Pruslin  [1],  in  a  simulation  of  the  receiver  noise 
effects  on  monopulse  performance,  has  observed  that  the  estimator  becomes 
biased  for  low  values  of  the  SNR.  Browne  [2],  in  an  interference -free  analysis, 
has  calculated  this  bias  exactly  for  all  values  of  the  SNR. 

Closely  related  to  the  bias  introduced  by  fading,  is  what  we  shall  refer 
to  as  the  azimuthal  bias  that  results  when  the  interferer  and  the  target  signal 
sources  are  located  at  different  azimuths.  Finally,  the  presence  of  inter¬ 
fering  signals  causes  the  variance  of  the  estimate  to  increase. 

In  this  note  we  shall  study  mainfceam  and  sidelobe  interference  and 
noise  effects  simultaneously  and  use  an  analysis  similar  to  Browne's  to 
compute  an  exact  expression  for  the  estimator  bias.  We  specialize  our  results 
to  real  antenna  patterns  and  obtain  an  expression  for  the  bias  that  holds  for 
the  high  and  low  SNR  situations,  hence,  we  obtain  the  azimuthal  bias  and  the 
fading  bias  simultaneously.  Numerical  results  are  given  for  some  typical 
operating  conditions  and  it  is  shown  that  the  bias  cm  be  a  large  fraction  of 
a  beamwidth. 

Some  further  analytical  results  are  obtained  for  the  case  of  sidelobe 
interference  for  real  antenna  patterns.  In  this  case  too,  it  is  possible  for  an 
ATCRBS  (Air  Traffic  Control  Radar  Beacon  System)  interferer  that  is  located 
in  a  large  near -in  sidelobe  to  cause  bias  errors  that  are  a  large  fraction  of 
a  beamwidth. 


In  addition  to  introducing  a  deterministic  bias  to  the  azimuth  estimate, 
interference  causes  the  random  error  to  have  an  increased  variance. 

We  obtain  a  general  expres  sion  for  the  variance  that  applies  to  the  mainbeam 
and  sidelobe  interference  cases  using  an  analysis  similar  to  that  developed 
by  Sharenson  [3]  who  analyzed  the  noise -only  case.  Numerical  results  are 
given  that  show  that  the  variance  also  depends  critically  on  the  relative 
phase  between  the  target  and  interference  signals. 

These  general  expressions  for  the  bias  and  the  variance  of  the  mono¬ 
pulse  estimate  provide  the  tools  needed  to  thoroughly  understand  the  effects 
of  interference.  These  are  essential  in  the  evaluation  of  signal  processing 
techniques  to  overcome  the  effects  of  interference.  The  leading  candidate 
at  the  moment,  is  the  monopulse  interference  detector  and  data  editing 
scheme  followed  by  an  outlier  test.  Although  this  is  the  subject  of  a  sepa¬ 
rate  paper  [4]  the  results  of  the  present  study  are  the  background  for  the 
succeeding  analysis. 

2.  INTERFERENCE  EFFECTS  ON  ESTIMATE  BIAS 

We  restrict  our  attention  to  sum-difference  (even- odd)  monopulse  pro¬ 
cessing  that  is  performed  on  a  sampled -data  oasts.  Assuming  mixer  pre¬ 
amplifiers  at  the  output  of  the  sum  and  difference  beams,  the  received  sig¬ 
nals  in  the  presence  of  interference  are  modelled  by 


i  *  1,  2 


Yi  =  V 


3<P. 


Gi(V  +  V 


3«Pt 


+  n. 


<D 


where  y^,  refer  to  the  outputs  of  the  sum  and  difference  beams;  Ag,  cpg, 

■» "  'N 

6g  and  Aj«  9j  are  the  amplitude,  phase,  and  azimuth  of  the  target  and  interference 
respectively;  Gj(  ),  G2(  1  are  the  antenna  patterns  of  die  sum  and  difference  beams, 
and  these  may  be  complex  in  general;  n^,  are  the  independent,  zero  mean  Gaussian 
noise  samples  due  to  the  mixer  preamps  whose  real  and  Imaginary  parts  have  vari¬ 
ance  a2. 

An  azimuth  estimator  that  is  often  used  in  practice  is  given  by 


the  relation 


E(8)  =  Re(y2/yi)  (2) 

where  the  monopulae  function,  E(9),  is  simply  the  normalized  difference 
pattern, 

E{6)  =  G2(6)/G1(9)  (3) 

In  most  cases  of  practical  interest,  this  function  is  well  approximated  by 
a  linear  characteristic,  hence  we  can  write 

E(9)  =  k(8/eB)  (4) 

where  9_,  refers  to  the  3  dB  beamwidth  of  the  antenna  patterns  and  k  is  a 
standard  parameter  that  arises  in  the  characterization  of  monopulae  systems. 
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Usually  1  -  k  -  2  with  k  =  1.  5  being  a  typical  value.  When  thi3  approximation 
is  valid,  the  azimuth  estimate  is  readily  generated  as 


(5) 


We  will  incorporate  6  into  our  definition  of  8  and  hence  express  all  our 

x3 

results  in  units  of  a  3  dB  beamwidth. 

Our  goal  is  to  compute  the  mean  and  variance  of  (5)  when  inter¬ 
ference  is  present  in  the  sum  and  difference  signals  as  given  by  (1).  We 
begin  by  letting 

jo  jcDj 

li.  =  A  G.(9  )e  +  ArG.(fi)e  =  V.e*  (6) 

x  s  s  III  1 


Since  m  are  complex  Gaussian  variates,  they  can  be  expressed  as 

=  N.  exp  ja!  where  are  Rayleigh  and  a!!  are  uniform  random  variables. 
Using  these  definitions  in  (1)  we  have 


y.  =  V.e 
i  i 


j°i  is. 


+  N.e  =  e 

i 


V.e  +  N.e' 

i  i 


(7) 
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V1  +'2v1N1co»ol  +N1 


(8; 


Following  Browne  [2]  we  first  average  (8)  over  ct^.  Since  this  is  uniform  on 
(0,  2tt),  we  have 


A  V2,cos(82  -  $  )  (Vj  +  NjCoso^)  +  V2N1sin(82  -  p^sinofj 

E  (k0)  =  — — —  - —  1  a  g - - - - 

®2  V1  +Ni  +  27  N  cosffj 


(9) 


Next  he  a<  erages  over 


In  the  Appendix  this  is  shown  to  give 


E 

al°2 


A 

(k  e ) 


0  if  Nj  >  Vj 

v2 

cos  02  -  pj)  if  Nj  <  Vj 


(10) 
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Since  Nj  is  the  magnitude  of  a  complex  Gaussian  random  variable  having 
2 

variance  la  ,  it  has  the  Rayleigh  distribution 


N 


1 


P(Nj)  =  —  exp 
rlc 


N. 


4  a 


for  N  =  0 


(11) 


Then  averaging  (10)  over  gives 


E(k6)  =  cos  ($2  - 


a  u 


I  »! 


■N. 


exp 


4a 


2  |dNl 


~  COS  (b2  -  Bj) 


•v2 


i  -  e*p  i  2 

2a 


(12) 


It  should  be  noted  that  when  interference  is  absent,  =  0 

and  from  (6)  we  see  that  =  g^  an<l  (12)  reduces  to  the  same  expression 

obtained  by  Browne.  However,  we  are  now  in  a  position  to  determine  the 

effects  of  interference  on  the  azimuth  bias.  To  do  this  we  need  to  evaluate 

Vj,  V2/V ^  and  8 ^  from  (6). 

Before  performirg  this  evaluation  we  first  note  that  the  target  of 

interest  lies  wi*.hin  the  mainbeam  of  the  antenna,  hence  we  may  assume  that 

G.(0  )  and  Go(0  )  are  real  functions  of  9  .  Whereas  the  sum  and  difference 
l  s  L  s  s 

beams  are  in  phase  over  the  antenna  beamwidth,  small  errors  in  the  ampli¬ 
tude  and  phase  tapers  render  them  complex  ir,  the  sidelobes.  Since  the 


6 


interference  can  be  located  in  the  mainbeam  or  within  the  sidelobes  of  the 
antenna,  then  G^(8j)  must  be  considered  to  be  complex  in  general.  To  make 
this  explicit  we  write 


G.(9j)  =  (Qj)  exp'f.ftj) 


(13) 


Then  substituting  this  into  (6)  we  can  easily  show  that 

V  |AsGi(0s)  +  2AsAIGi(9s)Ai(9I)CO8^+¥ei)^  +AK2(6I) 


1/2 


(14) 


-1  VVV  rtnfa+yej) 

e.  =  cs  +  tan  A  G.  (9  )  +  A  A  (9.)  cos  fcD+  t  (0T)1 

SIS  Jll  L  llJ 

where  cp  =  c.*^  -  co  is  the  relative  phase  between  the  target  and  interference 
signals.  Using  these  expressions  in  (12)  and  performing  straightforward 
but  tedious  trigonometric  manipulations  we  can  show  that  the  first  moment 
of  the  azimuth  estimate  is 


A  A 


[ 

.  2  2 

l 

A  G  (9  r 

1  -  exp 

8  I  s 

2 

( 

2(7  L 

A  (9) 

wp  TyiT  coa 


f2(9j) 


+  P 


Gi<V  J 


(16) 


where  p  =  Ar/A  represents  the  interference-to-signal  ratio  (ISR)  prior 

i  8 


to  any  processing  and 


ki  = 

o 


G-(0  )  ?  A  A 

_L^  +  ^_2._icos(T2-V  +  p| 


rG2<8s»  Al 


1 


A2  G2^V  A2  A1 

2  2’  2+2G(0)GG  cos^2~^l^t 

G,  (@  )  G,  G,  l' V  u)  U1  i 

1  s  1  1 


A  A2 

1  2  1 
1  +  2p  — —  cos  (cc»+  Y  )  +  p  r* 

1  Gx 


(17) 


where 


Y  =  tan 


-1 


G?(9S)  Al  sin  'Kj  +  G^)  A2  sin  *2 
C2  $g)  Ax  cos  ^  +  G^Sg)  A2cos  1!  ^  J 


(18) 


These  are  the  key  expressions  we  need  to  fully  understand  the  deterministic 
aspect  of  the  azimuth  error  introduced  by  interference. 

3 .  SPECIAL  CASE :  REAL,  ANTENNA  PATTERNS 

The  results  that  have  been  obtained  to  this  point  are  free  of  any 
approximating  assumptions  and  apply  to  the  general  case  of  complex  antenn« 
patterns.  Unfortunately  it  is  difficult  to  draw  ara^tical  conclusions  from 
these  expressions  because  of  the  large  number  of  parameters  that  would  have 
to  be  examined  simultaneously.  Results  for  the  general  case  will  be  obtained 
using  a  "omputer  simulation.  There  is  one  case  of  considerable  practical 
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importance  that  we  can  explor  e  analytically.  This  is  the  situation  in  which  the 
antenna  patterns  are  real,  which  is  a  reasonably  accurate  model  for  the  high 
near-in  sidelobes  of  the  antenna.  For  this  case  we  assume  that 


A.(0  )  sin  '{ (6  )  0  i  =  1 ,  2 

i  I  '-I 


(19) 


and  then  define 


G.  (9J  =  A.  (0  )  cos  ¥.  t  9  ) 
i  I  1  i  I 


(20) 


which  we  note  are  real  functions  for  all  values  of  9^..  Using  this  assumption, (16) 
and  (17)  become 


1  -  exp {  - 


2  2 

AsGl(9s) 

2/ 


,  G1(9I>  2G1<91» 

1  +  2p  . . .  cos  o+  p 


G  (9  ) 

1  s 


G2(0  ) 

1  s 


(21) 


wue  re  now 


A 

k  9  - 
o 


G^(0  ) 

-  ^  s 

G_(0  )  G.(0L)  GJfl.) 

2  s  i  i  2  1 

2G2(6I)  <W 

G  (0  ) 

1  s 

G  ( 0  )  G  ( fl  )  G  (9  ) 

u  1  s  Is  lSJ 

cos  cp  +  o  G  G  (0  ) 

IS  Is 

Gl(9I)  2G1<9I> 

1  +  2p  ~r  cos  rp  +  p  — r - 


G  (9  ) 

1  s 


Gl(9s> 


(22) 
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x 


We  recall  that  p=  A  /A  represents  the  in  te  r  f  cr  en  ce~  to  -  £ ignal  ratio  (ISR)  prior 

I  s 

to  any  processing.  The  measured  ISR  at  the  outputs  of  the  antenna  ports  is 
given  by 


AC.te)  0,(6.) 

vw'"w 


Although  both  definitions  of  ISR  are  important  parameters  of  the  system,  it 

is  convenient  to  rewrite  (21)  and  (22)  in  terms  of  p  .  In  this  case  the  mean 

o 

value  of  the  azimuth  estimate  is  given  bv 


9=  $  1  1  -  exp 


'As“Gl{Gs)  _  .2. 

5  (1  +  2  p  cosec+p  i 

£■  O  o 

2  cr 


where  now 


A 

k  e  = 

o 


G2<V 

r<yv  ‘vv  i 

2  G2(V 

G1<V  ‘ 

3 

o 

LG1(V  G1(9I)  J 

c08tr+0o  OjfSj) 

1  +  2p  cos  a,  +  p 
o  o 


It  is  interesting  that  the  performance  depends  only  on  the  monopulse  function 


E(9)  =  G2(8)/Gj(9).  In  most  cases  of  interest  the  target  will  be  located  within 


the  3  dB  beamwidth  of  the  antenna.  Therefore,  using  (3)  and  (4)  we  can  write 


where  6  is  measured  in  3  dB  beamwidths.  In  the  case  of  mainbeam  inter- 
8 

ference  this  linearity  property  will  also  be  valid,  but  in  more  general  cases 
E(9j)  will  o imply  take  on  the  values  of  the  normalized  difference  pattern.  To 
handle  this  situation  we  define  an  equivalent  azimuth 


#*w 

6, 


1  G2(V 

k  G^) 


(27) 


and  note  that  whenever  8^  is  within  the  3  dB  beamwidth  of  the  antenna  6^.  s  8^. 

In  Figure  1(a)  we  have  plotted  the  equivalent  azimuth  for  a  monopulse  antenna 
configuration  that  has  a  4°  beamwidth  and  achieves  -20  dB  peak  sidelobe  level 
on  both  the  "sum'1  and  "difference'1  beams.  Assuming  that  9j  is  uniformly 
distributed  in  (-t;,  tt),  Figure  1(b)  shows  the  probability  distribution  function 
of  the  equivalent  azimuth.  This  shews  that  "most  of  the  time"  the  equivalent 
azimuth  will  be  less  than  2  beamwidths,  hence,  for  the  purpose  of  analysis, 
we  can  study  the  cases  of  xnainbeam  and  sidelobe  interference  simultaneously. 
Then  using  (26)  and  (27)  in  (25)  we  can  express  the  first  moment  of  the  azimuth 
estimate  as 


A  A 

8  =  9 


1  -  exp 


2  2 

A  G  (0  ) 

— ~ — r — —  (l  +  2p  cosec+o  ) 
-  Z  o  o 


(28) 


a  9  +o  (6  +9T)cosm+0  9 
A  s  o  s  I  o  I 

6  =  - - 

o  ,  _  L 

1  +  2p  cos  cp+  o 
o  o 


(29) 
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PROBABILITY  DISTRIBUTION  EQUIVALENT  AZIMUTH 

FUNCTION  OF  EQUIVALENT  l  dig) 


EQUIVALENT  AZIMUTH  (  Beomwidths) 


Fig.  lb.  Probability  distribution  function  of  equivalent  azimuth. 


It  is  from  this  equation  that  we  can  begin  to  make  some  general 
statements  regarding  the  performance  of  the  azimuth  estimator.  We  shall 
first  show  how  the  results  are  consistent  with  those  previously  obtained  and 
then  go  on  to  consider  cases  that  have  not  yet  been  explored  in  the  literature, 
(a)  Case  1:  No  Interference 

When  there  is  no  interference,  AT  =  0,  hence  0  =0  and  from 

1  o 

(28)  and  (29)  the  average  value  of  the  estimate  is 

(30) 

This  is  the  result  obtained  by  Browne  [2]  and  describes  analytically  Pruslin's 

2  2  2 

observation  [1]  that  as  the  SNR,  Ag  Gj  (6  g)/2 o  .decreases,  the  monopulse 
estimate  is  biased  towards  the  antenna  boresight. 

We  have  tabulated  (30)  in  Table  1  from  which  it  can  be  seen  that 
for  SNR's  above  8  c(B,  the  bias  can  be  considered  to  be  negligible.  Since  in 
the  Air  Traffic  Control  (ATC)  system  the  SNR  will  be  at  least  20  dB,  this 
effect  will  oe  insignificant  unless  fading  occurs,  as  will  be  discussed  in  the 


next  section. 


Table  1.  Variation  of  Bias  with  SNR 


I - 

SNR  (dB) 

a>"  ^ 

•  CD 

l<<P 

0 

0.  37 

3 

0. 135 

3.6 

9.  2 

4.  76 

0.05 

6.6 

0.01 

7.25 

0.  005 

8.39 

0.001 

8.8 

0.  0005 

9.64 

0. 0001 

(b)  Case  2:  Fading 

If  the  target  and  interference  signals  arrive  from  the  same  azimuth 
(multipath  from  a  flat  earth  for  example)  then  6^.  =  =  6g  and  (23)  and  (29) 

reduce  to 
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(31) 


This  demonstrates  the  existence  of  a  fading  bias  that  can  occur  even  at  large 
SNR  when  the  interference  signal  occurs  out  of  phase  with  the  target  and 
tends  to  cancel  its  energy.  From  Table  1  we  conclude  that  the  bias  due  to 
fading  will  be  negligible  as  long  as  the  effective  SNR  is  greater  than  8  dB. 


(32) 
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For  a  20  dB  SNR,  this  shows  that  fading  could  become  a  problem  only  if 
.75<  P0  «  1*25- 

(c)  Case  3:  Azimuthal  Bias 

Next,  we  consider  the  case  in  which  the  SNR  is  large  and  the 

interference  is  of  low  enough  power  that  the  exponential  term  in  (28)  is 

negligible.  Then  (29)  reduces  to 

^  2 

T1  0  +  (8  +  0  )  o  cos  CD  +  p  0T 

0  =  -5 - S - L_° - J.  ,33 
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This  result  is  sufficiently  general  to  describe  the  cases  of  mainbeam  and 

~  A 

sidelobe  interference.  Notice  that  when  0  =  0T  =  0T  then  8=0  and  the 

s  I  I  s 

estimator  is  unbiased.  Therefore  a  bias  is  obtained  only  when  0  =£  6  T , 

si 

namely  when  the  target  and  interference  signals  are  separated  in  azimuth. 
It  is  for  this  reason  that  we  refer  to  this  effect  as  the  azimuthal  bias. 

To  obtain  some  physical  appreciation  for  the  behaviour  of  the 
estimator  we  consider  a  specific  case  of  interest.  We  assume  the 


target  is  located  on  bore  sight  at  20  dB  SNR.  Therefore,  0  =0  and 

s 

2  2  .2 

A  G  (0)/2a  =  100.  We  next  let  the  interferer  be  located  at  the  edge  of  the 

S  1 

'  >  dB  beamwidth  so  that  0^  =  8^  =  .  5.  To  within  a  good  approximation  the  sum 

2 

beam  is  quadratic  over  the  beamwidth,  hence  G^{0)  =  1-1.78  .  Since  the 
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output  interference -to-signal  ratio  is  p  =  pG  (8  )/G  (0  )  where  p  =  A  /A  , 

O  i  x  X  8  a  S 

we  sc  3  that  p  =  .  7C7p  Equation  (33)  is  plotted  as  a  function  of  the  relative 
o 

phase  angle  for  various  values  of  p.  The  results  are  shown  in  Figure  2. 
Equation  (28)  was  also  plotted  but  there  was  no  discernable  difference  in  the 
resulting  curves  at  least  for  the  values  of  p  wc  used,  hi  Figure  3  these 
equations  were  again  plotted  for  p^  nearer  to  unity  and  the  fading  effect  can 
be  observed  at  the  out  of  phase  condition.  Since  there  is  only  a  small  range 
of  values  for  which  this  effect  is  observable  and  since  its  effe:t  is,  if  anything, 
beneficial,  -we  shall  neglect  the  fading  bias  in  the  rest  of  our  work  and  restrict 
our  attention  to  the  azimuthal  error  as  described  by  (29). 

The  bias  in  the  estimate  can  be  written  as 


A  A  ~  p  +  cos  CD 

b(co)  i  0  -  8e  *  (0.  -  9  )  -2 - - - r-  (34) 

1  +  2p  cosc0  +  p 
*o  o 

where  pQ  and  0j  are  given  by  (23)  and  (27)  respectively.  Therefore  for  small 
values  of  pQ  ,  \  -  9  ,  while  for  large  values  §  =  0j  which  shows  how  the  inter- 
ferer  "captures"  the  estimate.  For  intermediate  values,  curves  like  those 
shown  in  Figures  2  and  3  are  obtained  which  demonstrate  the  so-called  scin¬ 
tillation  effect  which  is  a  term  used  to  describe  the  fact  that  the  azimuth 
estimate  lies  outsiie  the  azimuth  interval  (0g>  0j). 


D ASS  AZIMUTH  ESTIMATE  (Beomwidths) 


Fig.  2.  Azimuth  estimate  in  the  presence  of  interference. 
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(d)  Case  4:  Effects  of  Relative  Phase 

It  is  reasonable  to  model  the  phase  a3  a  uniformly  distributed 
random  variable  on  {0,  2rr).  Then  averaging  over  phase,  it  is  easy  to  show 
that 


2tt 

JL  f 

2rr  / 


b(cp)  dec  = 


if  p  <1 

rQ 


if  p  >1 
c 


(35) 


In  a  practical  situation,  the  key  issue  is  how  correlated  the  phase  is  from 
sample  to  sample  within  a  reply.  For  multipath,  the  correlation  time  can 
be  several  seconds  because  the  phase  relationship  depends  on  the  relative 
paJi  lengths  between  the  direct  and  multipath  signals  and  for  typical  aircraft 
speeds,  these  change  slowly.  For  ATCR3S  interference,  the  phase  difference 
will  be  independent  from  sample  to  sample  since  the  transmitter  tube  is  in¬ 
coherent  from  pulse  to  pulse.  However,  since  there  are  relatively  few  inter¬ 
ference  samples  in  any  one  reply,  the  inherent  averaging  due  to  phase 
cannot  reliably  be  exploited  and  we  must  face  the  possibility  of  having  to 
deal  with  the  large  errors  shown  in  Figure  2. 


(a)  Case  5:  Sidelobe  Interference 

Although  all  of  the  results  given  so  far  are  applicable  to  mainbeam 

and  sidelobe  interference,  there  are  some  remarks  that  are  worth  noting  for 

the  latter  case.  In  the  case  of  ATCRBS  sidelobe  interference,  it  is  possible 

that  Ar/A  can  be  much  larger  than  unity.  Let  us  suppose  for  example  that  the 
X  s 

DABS  target  is  at  maximum  range  (100  mi)  and  the  ATCRBS  interfe?’er  is  at 
the  same  power,  but  at  a  close  in  range  (10  mi  say).  Then  the  20  dB  decrease 
in  ISR  due  to  the  antenna  sidelobe s  is  compensated  by  the  20  dB  increast  in 
power  due  to  the  range  difference.  Then  oq  can  be  near  0  dBand  since  the 
equivalent  azimuth  is  of  the  order  of  1  or  2  beamwidths  then,  as  we  have 
seen,  large  bias  errors  can  be  expected.  Therefore  we  conclude  that  if  a 
strong  ATCRBS  signal  is  being  received  in  a  near-in  sidelobe,  the  monopulse 
azimuth  estimate  can  have  a  large  bias  that  exhibits  the  same  dependence  on 
the  relative  phase  between  the  DABS  and  ATCRBS  transponders  that  was 
shown  to  exist  for  mainbeam  interference.  Furthermore,  it  is  clear  that  very 
low  sidelobes  is  not  a  sufficient  means  of  eliminating  the  effects  of  this  inter¬ 
ference.  Therefore,  at  least  one  additional  level  of  data  processing  will  be 
needed  to  improve  the  quality  of  the  azimuth  estimate  for  these  cases.  This 
is  referred  to  as  monopuise  d'.ta  editing  and  will  be  discussed  in  a  later 
paper  [4], 
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For  sidelobe  muitipath,  on  the  other  hand,  we  can  reasonably  expect 
that  the  reflection  coefficient,  A^A^  ~  .  707.  For  antenna  patterns  with  20  dB 
peak  sidelobe  levels,  this  puts  less  than  .  07,  From  our  studies  so  far  we 
know  that  the  bias  error  peaks  when  the  direct  and  multipath  signals  are  it  out 
of  phase.  Using  (34)  the  bias  error  can  be  bounded  by 

Wp)|S  T^-lV  esl~  OV  V  (36) 

where  the  last  approximation  follows  from  the  fact  that  is  quite  a  bit  less 

than  unity.  For  a  DABS  target  on  bores ight,  6  =0  and  then  using  (23)  and  (27) 

8 

the  bound  becomes 


„  ,  A  G2(0t) 

ir  ~r  w 


(37) 


which  shows  that  the  bias  error  depends  on  the  sidelobe  level  of  the  differ¬ 
ence  pattern  relative  to  the  sum  beam  gain.  For  the  example  studied,  k=  1.  5, 
A. /A  -  .  707  and  G-, ( 0T)/  G,  (8  )  =  .  1 ,  which  yields  a  peak  bias  error  .  047 

Is  ills 

3  dB  beamwidths. 


4.  INTERFERENCE  EFFECTS  ON  ESTIMATE  VARIANCE 

In  this  section  we  will  compute  the  variance  in  the  monopulse  azimuth 
estimate.  Sharenson  [3]  has  studied  this  p’-oblem  for  the  case  when  the  inter¬ 
ference  consisted  only  of  receiver  noise.  We  shall  extend  his  analysis  to  in¬ 
clude  the  effects  of  mainbeam  and  sidelobe  interference.  We  begin  with  the 
equivalent  signal  model  formulated  in  (1),  (2),  and  15).  Tne  monopulse  estim¬ 


ator  of  interest  is 


(38) 


*  _  i  Re(yiy2> 

k  .  ,2 
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which  gives  the  azimuth  estimate  in  3  dB  beamwidths.  i  rom  (1)  and  (6), 


Y:  «  Uj  +  n. 


(39) 


so  that 


n 


1  ,  2 


h,r -m1 


(40) 


2,.  2 


where  the.  last  approximation  holds  provided  the  equivalent  SNR,  |(j  |  /2a 
is  large  enough,  typically  a  12  dB,  Then  the  moments  of  (38)  can  be  com¬ 
puted  by  evaluating  the  moments  of 


z  =  Re(y1y2  ) 


(41) 


Since  the  noise  terms  n^  and  are  independent,  the  first  moment  is  simply 


z  =  Rel 


(uiu2) 


(42) 


and  hence  the  first  moment  of  the  azimuth  estimate  is  given  approximately  by 


J  1V2  , 

8  =  -  —  cos  (0^  -  3,) 


(43) 
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which  is  the  same  result  we  obtained  in  the  last  section  when  the  effective 

SNR  was  greater  than  8  dB.  Proceeding  further,  the  second  moment  of  (4i] 

1  *  1 

can  be  calculated  using  the  identity  Re(x)Re(y)  =  —  Re(xy  )  +  — Re(xy.)  In  this 

$ 

case  we  let  x  =  y  =  y^y^  an<*  th®11 


2  I  ,  ,2  .  1 2  1  2*2 

z  =  2,yl!  ,y2f  +2Reyly2 


2  2  2 

From  (39)  and  using  the  fact  that  Jn.  (  =  2j  ,  n  =  0,  we  obtain 


ii2  it2  „  4 

|y.|  =  lu.f  + 


(45a) 


2  2 
yi  =  Ui 
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Combining  (42),  (44)  and  (45)  we  can  show  that  the  variance  of  z  is 


-2  / ,  2  .  2\  2  4 

=  ( luj!  +  lu2  1  |cr  J-  4cr 


Finally  we  use  (40),  (46)  and  (38)  to  show  that  the  variance  of  the  azimuth 


estimate  is 


a  1  a 
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Sines  liji  =  V  »  then  (47)  could  be  express. ?d  in  terms  of  the  complex  antenna 

pattern  parameters  by  using  (14).  Unfortunately  the  resulting  expressions 

are  complicated  and  in  order  to  obtain  some  meaningful  analytical  statements 

it  is  not  fruitful  to  retain  the  most  general  equations.  Therefore,  we  shall 

follow  the  direction  of  the  proceeding  section  and  restrict  our  attention  to  the 

case  of  real  antenna  patterns.  In  addition  we  note  that  (47)  describes  the 

estimate  variance  only  when  the  effective  SNR  is  large.  This  means  that 

.  ,2  2 

lUj  !  /^a  »1,  hence  we  can  also  neglect  the  effects  of  the  second  order 
SNR  term  that  appears  in  (47).  Under  these  conditions  we  can  then  use  (14) 
for  real  antenna  patterns  and  show  that  the  variance  reduces  to 
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where,  as  in  the  case  of  the  bias  effects,  0Q  and  9j  are  given  by  (231  and  (27). 

In  the  next  few  paragraphs  we  shall  consider  some  more  specialized 

cases. 

(a)  Case  1:  No  Interference 

When  there  is  no  interference,  Ar  =  0,  hence  p  =0  and  from  (48) 

I  o 


-.ve  see  that  the  variance  is  simply 


(49) 


which  is  the  result  obtained  Sharenson  [3]  and  others,  and  shows  that  the 
variance  decreases  with  SNR  and  the  slope  of  the  normalized  difference 
pattern  and  increases  with  the  degree  to  which  the  target  is  off  boresight. 

(b)  Case  2:  Fading 

If  the  target  and  interference  arrive  from  the  same  azimuth 
(multipath  from  a  flat  earth  for  example)  then  6^  -  6^  =  0g  and  (48)  becomes 


7  o 

A  2  1+k' V  . 

Var(0)  =  ♦  - j  (50) 

A  G  (9  )  k  l+2p  cos  cp+ p 

sis  o  o 

which  shows  that  when  multipath  fades  occur  there  can  be  a  reduction  in  the 
effective  SNR  which  in  turn  leads  to  an  increase  in  the  variance.  In  the  last 
section  we  found  that  a  bias  error  was  also  introduced  in  the  fading  situation, 
(c)  Case  3:  .Azimuthal  Variance 

When  the  target  and  interferer  are  at  different  azimuths  we  see 
from  (46)  that  the  variance  will  be  further  increased.  Typical  results  are 
plotted  in  Figure  4  for  the  case  of  a  target  on  boresight  at  20  dB  SNR 
and  an  interferer  located  at  the  edge  of  the  3  dB  beamwidth.  We  see  that  the 
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most  significant  degradation  in  performance  occurs  at  the  out-cf-phase  con¬ 
dition.  However,  the  magnitude  of  the  errors  is  much  smaller  than  the  con¬ 
tribution  due  to  the  bias,  and  furthermore  since  the  errors  are  random, 
their  effect  can  be  further  reduced  simply  by  averaging  many  of  the  estimates. 
Therefore,  we  conclude  that  the  bias  error  is  likely  to  be  the  more  troublesome 
problem  for  ATC  direction  finding. 

5.  CONCLUSIONS 

An  exact  expression  for  the  bias  of  a  two  beam  monopulse  azimuth 
estimate  has  been  derived  that  describes  the  degradation  in  performance  due 
to  receiver  noise  and  ATCRBS  and  multipath  interference.  At  low  SNR  and 
high  ISR  the  bias  tends  toward  boresight,  although  for  all  practical  purposes 
the  net  effect  of  this  bias  is  negligible.  More  important  is  the  azimuthal  bias 
that  describes  the  so-called  scintillation  of  the  azimuth  estimate  that  has 
been  observed  in  many  monopulse  tracking  problems.  Depending  on  the 
values  of  the  SNR,  ISR  and  the  relative  phase  between  the  DABS  and  inter¬ 
ference  signals,  the  bias  can  be  quite  significant  even  when  the  interfering 
signal  arrives  through  a  low  level  sidelobe. 

A  first  order  analysis  was  used  to  obtain  the  variance  in  the 
azimuth  estimate  when  interference  is  present  in  a  noisy  background. 

Although  the  results  indicate  that  the  random  errors  will  nol  ’'e  insignifi¬ 
cant,  the  deterministic  bias  error  will  be,  by  far,  the  more  dominant  effect. 


.  \ 
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For  ISR's  less  than  -"*0  dB,  the  effect  of  the  interference  is  negligible 
As  this  quantity  increases,  the  bias  and  variance  increase  and  errors  that  are 
of  the  order  of  a  beamwidth  can  be  obtained.  For  ISR's  between  _+  10  dB  there 
is  a  strong  dependence  on  the  relative  phase,  with  a  peak  error  occurring  at 
the  out-of-phase  condition.  As  the  ISR  is  further  increased,  this  peaking 
effect  subsides  until  the  interference  completely  captures  the  monopulse 
processor  which  at  this  point  would  track  the  interference  target. 

Although  the  results  are  applicable  to  analyzing  the  effects  of 
ATCRBS  interference  and  multipath,  a  distinction  between  the  two  phenomena 
should  be  noted.  For  multipath,  the  direct  and  indirect  signals  are  at  the 
same  frequency  and  coherent  in  the  sense  that  their  relative  phase  may  be 
constant  during  several  seconds  duration.  For  ATCRBS  interference, 
however,  the  transponders  are  incoherent  from  bit  to  bit  and  possibly  for 
samples  within  a  bit  since  there  may  be  carrier  frequency  offset.  Therefore, 
in  the  ATCRBS  case,  there  will  be  averaging  that  can  be  exploited  to  reduce 
the  overall  bias  error,  in  which  case  the  effect  of  the  increased  error 
variance  would  become  a  more  important  effect.  From  a  processing  point 
of  view,  this  could  be  overcome  by  using  more  samples  to  form  the  azimuth 
estimate. 

In  a  separate  study  [4j  the  performance  of  the  maximum  likelihood 
interference  detector  has  been  presented  in  detail.  The  next  step  is  to  com¬ 
bine  the  results  of  both  studies  to  evaluate  the  data  editing  concept.  This 
will  determine  whether  there  is  any  promise  to  the  idea  of  introducing  an 
additional  level  of  data  processing  to  improve  the  overall  quality  of  the 


azimuth  estimate. 
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APPENDIX 


To  compute  the  average  of  (9)  with  respect  to  o'j*  a  uniformly 
distributed  random  variable  on  (0,  2tt)  we  have 


<k0)  =  cos(8?  "  B.) 
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Since  the  integrand  of  the  second  term  is  an  odd  function,  its  integral  is 
zero.  Using  standard  integral  tables,  Browne  showed  that  the  integral  in  the 
first  term  reduces  to 
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where  the  positive  square  root  is  understood.  Therefore 


0  if  N-  <  V 

V2 

~  cos{c^  -  pj)  if  Nj  >  Vj 

as  given  in  (10). 
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